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Abstract 

We present a computation of nucleon mass corrections to nucleon structure 
functions for polarized deep-inelastic scattering. We perform a fit to existing 
data including mass corrections at first order in m 2 /Q 2 and we study the effect 
of these corrections on physically interesting quantities. We conclude that mass 
corrections are generally small, and compatible with current estimates of higher 
twist uncertainties, when available. 
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1. Introduction 



Experimental information on deep-inelastic scattering of polarized leptons off dif- 
ferent kinds of polarized nucleon targets has become more and more accurate in the 
past few years |J-||. This accumulated knowledge, combined with recent theoretical 
progress in the computation of perturbative QCD quantities relevant to polarized 



deep-inelastic scattering [10], has allowed a next-to-leading order determination of 
the polarized parton densities, using data on the structure function gi(x, Q 2 ) which 
determines the cross section asymmetries in the case of longitudinally polarized lep- 
tons and nucleons in the Bjorken limit Reasonably good determinations of 
the strong coupling constant a s and of the axial vector coupling g& have also been 
performed |TJ] . 

A large part of experimental data in polarized deep-inelastic scattering are taken 
at relatively low values of Q 2 . In particular, Q 2 is usually around 1 GeV 2 for data 
points in the small-x region, which is particularly interesting because there the effects 
of Q 2 evolution are more important (data at even lower values of Q 2 are also available, 
but they are usually not included in perturbative analyses). In this kinematical re- 
gion, contributions suppressed by inverse powers of Q 2 can become important. These 
contributions can be of two different origins. There are power-suppressed terms aris- 
ing from the operator product expansion of the hadronic tensor W^ u . These terms 
originate from matrix elements of operators of non-leading twist (they are usually 
referred to as dynamical higher twists). Their effect is not controlled by perturbation 
theory, and it is very difficult to assess their importance. A second class of power- 
suppressed contributions originates from taking into account the finite value of the 
nucleon mass m in the kinematics of the leading-twist cross section. These correc- 
tions can be computed exactly, and have been studied in detail in the past |15| in 
the case of unpolarized deep-inelastic scattering. Knowledge of kinematic effects is 
of course necessary in order to extract information on dynamical higher twists from 
experimental data. 

The problem of calculating kinematic higher twist terms in polarized deep-inelastic 
scattering has been considered in ref. [16|.Q There, the reduced matrix elements a n , d n 
of the relevant operators in the OPE were expressed in terms of polarized structure 
functions, taking mass corrections into account; these expressions reduce to moments 

1 The results of ref. jl6| have been used in ref. pjj to compute target mass corrections to the 
Bjorken sum rule. 
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of the structure functions in the massless limit, but do not have a simple parton model 
interpretation in the case m^O. For this reason, the result of ref. ]TB| is not directly 
applicable in a full analysis of polarized deep-inelastic scattering data. 

In this paper, we calculate target mass corrections in polarized deep-inelastic 



scattering extending the analogous work of ref. |nj for the unpolarized case. This 
has the advantage of yielding the final result in a form which is appropriate for 
phenomenological applications, that is, we will express moments of polarized structure 
functions as functions of the reduced operator matrix elements. The details of our 
calculations are presented in Sect. 2. An important point is of course the interplay 
between dynamical and kinematical higher twists; in Sect. 2 we discuss this point in 
some detail. In Sect. 3, we use mass-corrected formulae in an analysis of existing 
data in the framework of QCD at next-to-leading order, and compare the results 
with those obtained in the massless limit; in particular, we will compare our results 



with those of ref. [II , where an estimate is given of the uncertainties coming from 
higher twist effects on the determination of physically interesting quantities. We will 
see that the effect of mass corrections is indeed within the uncertainties estimated in 
ref. \TM. Finally, we present our conclusions in Sect. 4. 



2. Calculation 

Our calculation follows closely the analogous one in the case of unpolarized deep- 
inelastic scattering, performed by Georgi and Politzer back in 1976 |15[ 



We will 



adopt the usual notation for the definition of kinematical quantities in deep-inelastic 
scattering, used for example in ref. [pl|]. In the formalism of the operator product 



expansion, the antisymmetric part T^ v of the forward amplitude relevant for deep- 
inelastic scattering, 



T^ip, q,a)=ij d 4 xe illx <p } s\T [J^x)J u (0)} \p,s> (2.1) 



is given by [[il| 




E 



<T/il.../i„_l 



p,s> 
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n — 1 



Act/ii —Hn-2 



n 



p,s> 



where J M is the electromagnetic current, and Q^ u Xa is defined as 
The operators j^ 1 - ^™- 1 are given by 



(2.2) 



(2.3) 



CjUl—A'n-l ^n — 1 



rt l,1/j 

K 



i = 1 



5 * * * ? 



F 01 D>* 1 ...D^- 2 F^- 1 



(2.4) 

(2.5) 
(2.6) 



Here _D M is the QCD covariant derivative, and the symbol [. . .]s means complete 
symmetrization in the indices a, fa, . . . , /i n _i; i 7 ^ is the usual QCD gluon tensor. 
The twist-3 operators are gi ven \yy 

-,Act/ii.../i„_2 -n— 1 



2,; 



pAo-/ii.../i n _2 -n— 1 



^ 5l x D CT D^...D^- 2 ^i) 



, i = l,..., 8; 



ipj 5 -f X D a D^...D^- 2 i) 



e aa ^F^D^...D^- 2 F> 



s> ' 



(2.7) 

(2.8) 
(2.9) 



where [. . indicates antisymmetrization with respect to A and a and symmetrization 
with respect to other indices. The coefficient functions E™ i (Q 2 ,a s ) and E2 i (Q 2 ,a s ) 
have been computed up to order a 2 in perturbative QCD The expansion of the 
forward scattering amplitude in powers of m 2 /Q 2 is independent of the perturbative 
expansion in powers of a s ; we can therefore perform our calculation at leading order in 
a s , and then insert perturbative corrections to coefficient functions in the final result. 
For a s = we have simply E lti = E 2 a — 1 for i = 1, . . . , 8, tj) and E^ G = E 2 ,g — 0. 



The matrix elements of the operators ( |2.4| - |2.5| ) can be written as 

<p,s\ R^-^- 1 \p, 8 >= -2ma n Ml^-^- x 



(2.10) 



(we have omitted the index i, which is no longer necessary at leading order). The 
tensor M^ 1 '"^ 1 is the most general rank-n symmetric tensor which can be formed 
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with one spin four- vector s and n — 1 momentum four- vectors p; furthermore, it must 
satisfy the tracelessness conditions 

giHlh MP"** = Q (2.11) 

for all pairs With these requirements, we find 

71-1 

= ( ~f [sp.. .p] s (m 2 Y (2.12) 

J ~° 3 

up to an overall normalization, which can be absorbed in the definition of the reduced 
matrix elements a n . The symbol 

g_-_g (2.13) 

3 

represents a product of j metric tensors g w ^ fe , with indices chosen among fix, . . . , \i n in 
all possible ways; the remaining n — 2j indices of Mf 1 "^" are carried by the symmetric 
product [sp . . .p] s . When the nucleon mass is neglected, only the first term of the 
sum, [s^p^.-.p^-^g, is retained; this is the standard result, used for example in 



ref. [IS 



Consider now the twist-3 operators of eqs. (|2.7| - |2.9| ). Their matrix elements can 
be written as 

<p,s\ R^-^~ 2 \p,s>= mdnM^ 1 -^- 2 , (2.14) 

where the tensor M2 must be antisymmetric in (A, a), symmetric in all other indices, 
and traceless. It is easy to prove that 



^-A(T/ii...At„_i _ n + 1 /„crj-|-AAl-Atn-2 „AjJ(7Ml-Atn-2^ | U \ (rfF j^-pti.../*»_a\ 

n n ^ 

(2.15) 

where 

XL 

2 (—1)3 (n _ i)\ 

n" 1 ^^ J , 9---9 P-..p(m 2 y (2.16) 

3=0 "' ^3^ 

is the most general rank-n symmetric, traceless tensor that can be formed with the 
momentum p alone. For m 2 = 0, one recovers the usual result 

<p,s\ Btf* 11 -^- 2 \p,s>= md n (s> A - sV) P 1 * 1 ■ --P^- 2 - (2.17) 

The reduced matrix elements a n , d n contain all the information on the proton spin 
structure; they are related to moments of polarized structure functions. Our next step 
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consists in obtaining these relationships in the general case m ^ 0. To do this, we 
compute explicitly the amplitude T^ u using our results, eq. ( |2.10| ) and eq. ( |2.14j) . We 
decompose T^ v into a twist- 2 and a twist-3 component, 

T£, = iP+Tj$. (2.18) 

Let us first consider the contribution of twist- 2 operators, 

2#> = 2ime uvXa q x E (^-] q m . . . q^ Mf^ 1 a n . (2.19) 

n odd J 

Using eq.( |2.12|) and recalling that x = Q 2 /(2p ■ q) we find 



2™ 2 \ 3 



x m 



p-q 



n odd 71 j=0 



Q 2 J j\(n-i- 2 j)\ 



s a + ( n -2j- if—V* 
(2.20) 

Following ref. |15|, we change summation index from n to I, with n = 21 + 2j + 1, and 
exchange the summation order of I and j. This gives 

(m 2 \ (2/+J + 1)! 



2im 



= — W E *- {2l+l) [ s° + k—v° ) E 



p-q 



1=0 



P-Q Jf^(2l + 2j + iy VQ 2 



We now define functions F a4 (x) by 

a n = f dy y n F a (y); d n = f dy y n F d (y). 
Jo Jo 



It is easy to prove that 



On _ 

n ./(i 

dn _ 
n Jo 



dyy n - v G a (y); 



n* jo 



fdy y^HM 
Jo 



1 dy y"- 1 G d (y); % = t dy y^ 1 H d (y), 



Jo 



my- 

(2.21) 



(2.22) 

(2.23) 
(2.24) 



where 



G a4 {x) = f 1 dyF a4 (y); H a4 {x) = C —G ad (y). (2.25) 
Jx Jx y 



These definitions allow us to perform the summation over j in eq. (|2.21|) ; in fact, we 
can write 

2im 



rp(a) 

fll/ 



p-q 
1 



e,,A CT g A f>- (2m) (2/ + l) (s° + 2l^-p° 
1=0 \ P ' 1 j 







dyy "MT,{-qt) m + 1){ ■ 



^2.26) 
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The sum over j can now be performed using the identity 



1 d n 1 



1 - z) n+1 n\dz n l- z n 



3=0 



[n + j)\ 



^2.27) 



with the result 
2im 



;i - nf-^-pA f dy - ^ t H a (y) 
p-g /jo m »v ^ 



(2.28) 

where we have defined n — 21 + 1. The same procedure applied to the twist-3 term 
yields 



T (d) = ^ £ 
M1/ p-q 



oo „i y™" 1 

E ^ n / ^ J y^fiV 
n=3,5,... J0 1 — rW 



(n-1) 



G d (y) s °-nH d {y)^p° 
p-q 



+ 2n y 2 m 2 G d (y) + H d (y) ^ 



Q 



2 -r y 2 m 2 

1 



^2.29) 



It is now possible to obtain the n th moments of the polarized structure functions 
gi and g^. The antisymmetric part of the hadronic tensor W^ v is defined by 



iW, 



-- -ImT A . 



[2.30) 



From the analytic structure of T^ u in the complex v = p- q plane, one can prove that 



9 oo i oo i 

T i = ~ E W dyy^lwT* = E x~ n \ dyy^ 1 2iW*. (2.31) 

71 71=1,3,... J0 71=1,3,... J0 



In other words, the coefficient of x n in gives twice the n moment of the hadronic 
tensor iW„„. 



On the other hand, iW^ u is usually parametrized in the following way: 



iW. 



im 



p-q 



gi(x,Q 2 )s° + g 2 (x,Q 2 ) s CT - 



1- s ^ 



p-q 



(2.32) 



We can therefore identify the n th moment of g\ + g^ and gi as twice the coefficients 
of s a and — p a (s ■ q)/(p ■ q) in eqs. (|2.28| , |2.29"1) respectively, thus obtaining 



<(Q 2 ) = / dy — 
Jo (i 



,71-1 



y 2 m 2\n+l 



n 2 H a {y) + 



2y 2 m 2 

~Q 2 



nG d {y) + n 2 H d {y) 
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(2.33) 



g n 2 (Q 2 



n[n 



,n-l 



dy 
o fl 



H d (y)-H a (y) 



y m 



Q 2 



H d {y) 



[2.U) 



Equations ( |2.33| , |2~34"|) are our main result. They express the moments of the polarized 
structure functions gi, g 2 in terms of matrix elements of the operators appearing in the 
light-cone expansion of the forward scattering amplitude, at all orders in m 2 /Q 2 . Ob- 
serve that when the twist-3 operator matrix elements d n are neglected, eqs. ( |2.33| , 12.34| ) 



obey the so-called Wandzura-Wilczek relation [20 

1 



n 



n 



g?(Q 2 ) (for d n = 0). 



^2.35) 



The familiar m = result [18] is easily recoverd by using eqs. (2.22-2~2i): 



m 



9i(Q ) = a n + I — 



92 (Q 2 



n — 1 



77 



m 



(dn ~ On) + O . 



(2.36) 
(2.37) 



The inverse Mellin transforms of eqs. ( |2.33|j2.3^ ) can be computed as in ref. 

d 2 d 



gi(x,Q x 



X dx 2+X dx 



2xm 2 d 
Q 2 dx 



-G d (0 
r 



92(x,Q' 



x 



dx 2 



x 



r£ I x 



H d (0 - H a (0 



where 



2x 



1 + 



Ax 2 m? 



(2.38) 
(2.39) 

(2.40) 



1+r' V" ' Q 2 ' 

An excellent test of the correctness of our calculation is a comparison with the results 
of ref. [r6|. We have checked that eliminating g\ and gi from eqs. (18,19) of ref. [16) 
using our eqs. ( |2.38| , |2.39"| ) (a non-trivial task) two identities are obtained. 

As in the unpolarized case, there is a well-known difficulty in eqs. ( p. 381 , ^39]) at 
x = 1: in fact, when x = 1 the structure functions should vanish for kinematical 
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reasons, while the RHS of eqs. ( |2.38|j2l39|) are clearly nonzero, since £(x — 1) < 1 



This problem was discussed in ref. [^TJ for the unpolarized case; the conclusion reached 
there is that in the large x region dynamical higher twist corrections become important 
and cannot be neglected any more. This is because the twist 2 + 2k contribution to 
the n th moment of a generic structure function has the form 

BUQ 2 )(^f)\ (2-41) 

where A is a mass scale of the order of a few hundreds MeV, and the coefficients 
Bkn{Q 2 ) have no power dependence on n, k or Q 2 . The crucial feature of eq. ( |2.41 



is the presence of a factor n k , which arises because there are at least n twist 2 + Ik 
operators of a given dimension for each leading twist operator of the same dimension. 
One can prove that the behaviour of structure functions in the x ~ 1 region is governed 
by moments n = (D(Q 2 /m 2 ); in fact, when x = 1 and m 2 /Q 2 << 1 we have 

2 

m 

f = 1 - 0* P.42) 

on the other hand, if we assume a (1 — £) a behaviour for the structure function, with 
a of order 1, its n th moment receives the dominant contribution from the region 

£ ~ ~ i _ °L. (2.43) 

a + n n 

Comparing eqs. ( |2.42| ) and ( |2.43| ), we obtain that the relevant moments for the x = 1 
region are of order 

aQ 2 

n = ^r 2.44 
as announced. Inserting this in eq. ( |2.41| ), one immediately realizes that the contri- 



bution of dynamical higher twists is no longer suppressed by inverse powers of Q 2 
when x is close to 1, and we cannot expect our result, eqs. ( |2.33| , |2.34"D , to hold in this 



region. A solution to this problem |22], |23| is that of expanding the result in powers 
of m 2 /Q 2 up to any finite order. In this way, the dangerous contribution of terms 
with large powers of m 2 /Q 2 is not included. The expansion remains reliable even 
when Q 2 is as low as m 2 , provided x is not too large; in fact, powers of m 2 /Q 2 always 
appear multiplied by an equal power of x 2 . The expanded result of course cannot be 
reliable at x ~ 1, but this would not be the case even without expanding in m 2 /Q 2 , 



since we are not including the contributions of eq. (|2.41| ), which are important in 
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this region. Therefore, we will perform our phenomenological analysis expanding our 
results, eqs. ( p.33| , |2.34] ), to first order in m 2 /Q 2 . We have 

m 2 n(n + 1) f __ AJ s ^ ( m A \ 

(2.45) 



9^) = «n + 



Q 2 {n + 2) 



o K na n+2 + 4d n+2 ) + O I — 



m 2 n(n — 1) r , , ... . /m 4 \ 



#2 (Q ) = rf„-«n + 

n 



(2.46) 



We can now use eqs. ( |2.36| , |2.37| ) to eliminate the matrix elements a n , d n from eqs. ( p. 45 



[2.46D in favour of the moments of g\ and g 2 at zero nucleon mass, which we denote 
by g™ and g% , respectively. We obtain 



g n dQ 2 ) = aUQ 2 ) + 



92(Q 2 ) = 9UQ 2 ) + 



n /^2N , m2 U ( n + !) 



Q 2 (n + 2) 2 

m 2 n{n — 1) 
Q 2 (n + 2) 2 



n + 2 



» + 4) ^V 2 (Q 2 )+4^ <? 2 "o +2 (Q 2 ^ 



(2.47) 



n + 2 



n- 



j92o 2 (Q 2 ) - 9io 2 (Q 2S 



n + 



° SO- (24f 



3. Phenomenology 



Analyses of polarized deep-inelastic scattering data in the context of QCD at next- 
to-leading order have been performed by different groups |IT|- ||14{| in the zero- mass 
approximation. In this section we will repeat the same analysis taking mass correc- 
tions into account, in order to establish their practical importance. The quantities 
which are directly measured in polarized deep-inelastic scattering experiments are the 
asymmetries 



A 



A,, 



a 



a 



0-T 1 1 1 



a 



(3.1) 



(3.2) 



a U + <jTt ' 

where the arrows refer to the orientation of the lepton and the proton spin vectors 
with respect to the beam axis. One can show that 



A ± = d{A 2 - (A,) 
A, = D{A 1 + V A 2 ) 



(3.3) 
(3.4) 
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where A\ and A 2 are virtual photon cross section asymmetries, and the coefficients 
D, d, r), C are fixed by the kinematics of the process; in the target rest frame 



E{l + eR)' 1 E-eE r Vl + £ 2e v ; 

where £7 (E') is the energy of the incoming (outgoing) lepton, 

and (9 is the lepton scattering angle. The ratio R is the usual quantity defined in 
unpolarized deep-inelastic scattering, namely 

F 2 ( 4x 2 m 2 \ , s 

The asymmetries A\ and A 2 are directly related to polarized structure functions 
through 

A 1 = (3.8) 

M ~ v^^r- (3 - 9) 

Therefore, both g\ and g>2 can be expressed in terms of the measured asymmetries 
j4x,v4|| and of unpolarized structure functions, through eqs. ( |3.3j , |3.4| ) and ( |3.8| . |3.9| ). 
When the nucleon mass is neglected, one has simply 

^ = A 1 = ^ = ^2x(l + R) (form = 0) (3.10) 

while in general this relationship involves both An and A±. 



The usual procedure for analysing data, adopted for example in ref. |Tj], must 
be modified in different aspects. First, the term proportional to m 2 /Q 2 in eq. (|3.7|) , 
which relates F\ to R and F 2 , must now be included. Secondly, one must perform 
a global fit of the measured asymmetries using eq. ( p.47|) (if the analysis is done 
in moment space), where the moments <?" are given in terms of moments of the 
coefficient functions and of the polarized parton distributions, which by definition 
are proportional to the matrix elements a n . A difficulty immediatly arises, because 
moments of the structure function g 20 also appear in eq. (|2.47|) ; therefore, it is not 
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possible to treat g\ and g 2 independently, as in the m = case. One could in principle 
circumvent this problem using experimental information on g 2 ; unfortunately, g 2 data 
available up to now are restricted to a very limited range of x and Q 2 , and are affected 
by large uncertainties [23|. We prefer here to follow a different strategy. We will 
perform fits to data in two different ways, characterized by two different assumptions 
on g 20 : either we simply set g 20 = 0, or we relate g 2 o to g w by the Wandzura-Wilczek 
relation, 

g^(Q 2 ) = -^g"o +2 (Q 2 )- (3.n) 

Notice that none of the two assumptions is theoretically justified: there is of course 
no reason to assume that g 2 Q vanishes, nor that twist-3 operators give a negligible 
contribution. However, both assumptions are consistent with presently available in- 
formation on g 2 . We will check that this procedure actually allows to make a reliable 
estimate of the effect of mass corrections. Notice that some of the experimental 
collaborations present values of the asymmetry Ai, while others give values of the 
combination of A\\ and A± which corresponds to g\jF\. The two quantities coincide 
for m = 0, as already observed, but they do not when mass corrections are included; 
so in the first case the asymmetry must be fitted by [gi — (4m 2 x 2 /Q 2 )g 2 ]/Fi, in the 
second case simply by g\/F\. 

The structure function g w is related to the polarized quark and gluon distributions 

by 

9io(x, Q 2 ) = [C NS ® Aq NS + Cs ® AS + 2n f C g ® Ag] , (3.12) 

where (e 2 ) = nj 1 Y^=i e fi ® denotes convolution with respect to x, and the nonsinglet 
and singlet quark distributions are defined as 

A( ins = E (rk ~ 1 + A ^ i) ' AS = i + A ^' (3 - 13) 
i=i \ \ e / j i=i 

where Aqi and Aqi are the quark and antiquark distributions of flavor i and Ag is 
the polarized gluon distribution. 



A first fit (called fit A in ref. fI3}| ) is performed with the initial parton distributions 
parametrized at Qq — 1 GeV 2 according to the conventional form 

A}\x,Ql) =Af fVf x a f(l-xff(l + V x 5 f) (3.14) 

where Af denotes Aq^s, AS or Ag; the factor Mf is chosen so that the first moment 
of Af is equal to rjf. We have fixed 

6e = 6 b = 1, 5 ns = 0.75, p g = 4, 7s = 7g (fit A). (3.15) 
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Fit A 




solid: m=0 




, dashed: m=0.93B GeV, g a0 =0 




/ \sdottcd: m=0.93B GoV, g z0 =g z ™ - 




*\ Q 2 =10 Gev s — 











Fit B 

solid: m=0 

dashed: m=0.93B GeV, g z0 = 
dotted: m=0.938 GeV, g zo =g Z o 



Figure 1: The structure function g\ for the proton at Q 2 = 1, 10 GeV 2 in 
fits A and B, with and without target mass corrections (parameters fixed at 
the m = values). 



In a second class of fits (fit B) the rise at small x is at most logarithmic: 
AS = 7V E r/ s x as (log 



NnsVns (log 1/ x) aNS + -f NS x (log 1/ x) 



(fit B), (3.16) 



Aq NS = 

Ag = Af g T] g [(logl/x)" 9 + 7 9 x (logl/x) 
The motivations for the choice of these particular parametrizations are discussed in 



ref. 14 



We have first performed fits A and B with fixed a s (mz) = 0.118 for m — 0, 
and with the obtained values of the fit parameters we have recomputed the structure 
function g\ using the mass-corrected formula, eq. (|2.47 ). The effect of target mass 
corrections turns out to be very small. As an example, in fig. [I] we show gi(x,Q 2 ) 
for the proton at Q 2 = 1 GeV 2 and Q 2 = 10 GeV 2 and m = for both fits A and 
B; in the same plot, we also show gi obtained including mass corrections, using the 
same values of the fit parameters found for m — 0, for both assumptions #20 = 
and #20 — 92o W - Differences among the three curves for Q 2 = 1 GeV 2 are sizeable 
only in the large- a; region, while at Q 2 = 10 GeV 2 the three curves are practically 
undistinguishable on this scale. 

We have then repeated the same fits for m = 0.938 GeV with both assumptions on 
920, #20 = and (720 = 92o W - The results for g\ (proton) are shown in fig. ||. Observe 
that also in this case the difference between the solid curves, which correspond to the 
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Fit A 

solid: m=0 

dashed: m=0.938 GeV, g 20 =0 
^\dotted: m=0.938 GeV, gz^gza* - 






\. Q a =10 Gev 2 — 








-3 


10~ 3 


10 _1 10 





Fit B 




solid: m=0 




dashed: m=0.938 GeV, g 31) = ~ 




dotted: m=0.938 GeV, gao = Sao W 




'^X ^-^---^=10 Gev a 




^^O a : 1 Gev 2 ^^>- 







Figure 2: The structure function g\ for the proton at Q 2 = 1, 10 GeV 2 in 
fits A and B, with and without target mass corrections (fitted parameters). 



m = case, and the dashed and dotted curves, corresponding to m / 0, are quite 
large at Q 2 = 1 GeV 2 for x > 0.2; when x is very large (x > 0.7), however, the effect 
is not physical because in this region the approximation under which mass corections 
have been computed is not reliable. However, these differences become negligible at 
higher vales of Q 2 because of the m 2 /Q 2 suppression factor, as one can see by looking 
at the Q 2 = 10 GeV 2 curves; this guarantees that the uncertainty of our calculation 
connected with the expansion in m 2 /Q 2 at first order has little effect on the results 
of our fits to data, since data at large x are taken at high values of Q 2 . The effect at 
small x, x < 0.01, is also rather large; however in this region there are few or no data 
and thus gio itself is affected by a substantial uncertainty. 

The results of fits A and B for m ^ 00 for some of the fitted parameters are 
shown in the second and third columns of table p]. The errors quoted in the table are 
statistical errors originating from experimental uncertainties only. 



In ref. |14j the uncertainty on from higher twist corrections was estimated to be 
±0.03. We can see from table [l] that target mass corrections actually induce deviations 
of this order from the m = value for g^. Deviations are always in the direction 
of making gA larger than in the m = case. We observe also that the values of g& 



2 The results in the first column of table [j] are slightly different from those of ref. 14 for two 
reasons: first, we are fixing (3 g = 4 in fit A, second, we are using more recent data sets for proton 
SMC data and neutron E154 data. 
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FIT A 




171 = 


920 = 


920 — 920 


9a 


1.167 ±0.045 


1.192 ±0.040 


1.200 ±0.043 




0.426 ±0.037 


0.408 ±0.027 


0.416 ±0.031 


% 


0.98 ±0.25 


0.81 ±0.35 


0.83 ±0.32 


a (10 GeV 2 ) 


0.19 ±0.04 


0.21 ±0.08 


0.21 ±0.07 


d.o.f. 


114-10 


114-10 


114-10 


x 2 


91.2 


86.9 


89.5 


xVd.o.f. 


0.88 


0.84 


0.86 


FIT B 




m = 


920 = 


920 — 920 


9a 


1.253 ±0.057 


1.292 ±0.056 


1.277 ±0.058 




0.455 ± 0.038 


0.423 ± 0.034 


0.428 ±0.037 


% 


1.40 ±0.32 


1.00 ±0.33 


0.99 ±0.35 


a (10 GeV 2 ) 


0.13 ±0.05 


0.18 ±0.06 


0.19 ±0.06 


d.o.f. 


114-11 


114-11 


114-11 


x 2 


86.8 


86.3 


90.1 


xVd.o.f. 


0.84 


0.84 


0.87 



Table 1: Results of fits A and B with fixed a s (mz) = 0.118. 
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FIT A 




m = 


#20 = 


920 — 920 


a s (m z ) 


0.118 ±0.005 


0.117 ±0.004 


0.120 ±0.003 




0.433 ± 0.039 


0.415 ±0.027 


0.423 ±0.028 


Vg 


1.04 ±0.45 


0.91 ±0.29 


0.85 ±0.40 


a (10 GeV 2 ) 


0.18 ±0.08 


0.20 ±0.09 


0.20 ±0.09 


d.o.f. 


114-10 


114-10 


114-10 


x 2 


94.4 


89.1 


90.9 


X 2 /d.o.f. 


0.91 


0.86 


0.87 


FIT B 




m = 


920 = 


920 — 920 


a s (m z ) 


0.123 ±0.003 


0.118 ±0.005 


0.121 ±0.004 


V^ 


0.448 ± 0.036 


0.407 ±0.036 


0.418 ±0.033 


Vg 


1.01 ±0.32 


0.73 ±0.33 


0.72 ±0.31 


a (10 GeV 2 ) 


0.14 ±0.07 


0.22 ±0.07 


0.21 ±0.06 


d.o.f. 


114-11 


114-11 


114-11 


x 2 


84.8 


86.3 


89.2 


xVd-o.f. 


0.82 


0.84 


0.87 



Table 2: Results of fits A and B with fixed gA = 1.2573. 

obtained with the two different assumptions for #20 are quite close to each other, thus 
suggesting that the assumed form of g>2o has a small impact on the final results. 

We have also performed fits A and B using a s as a free parameter, with fixed 
to its measured value = 1.2573. The results are shown in table Also in this 
case, mass corrections induce variations on a s (mz) which are compatible with the 
higher twist uncertainty of ±0.004 estimated in ref. [jT4f| . 

The values of the first moment of the quark singlet distribution, rfe, and of the first 
moment of the gluon distribution at Q 2 = 1 GeV 2 , rj g , are also shown in tables [l]-[| 
Also for these quantities, the assumption on #20 has little effect. The introduction of 
mass terms tend to give smaller values for and 7] g , but deviations from the values 
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Fit A 

solid: m=0 

dashed: m=0.938 GeV, ga = 
dotted: m=0.938 GeV, g 30 =E30 (n 





Figure 3: The quark singlet distribution at Q 2 = 10 GeV 2 in fits A and B, 
with and without target mass corrections. 



of the m = case are all within statistical errors. 

In tables we also present values of the singlet axial charge 



a (Q 2 



dx 



AE(x,Q 2 )-n 



2tt 



-Ag(x,Q 2 



(3.17) 



for Q 2 = 10 GeV 2 . Values of for m ^ are slightly larger than in the massless 
fits. The only case in which a sizeable devitation from the m = result is observed is 
fit B with fixed. Also in this case, however, the difference is compatible with the 
statistical error. 

Finally, we compare the polarized parton distribution functions for quark singlet, 
quark non-singlet and gluon obtained with and without mass corrections. They are 
shown, for Q 2 = 10 GeV 2 in figs. §, |]and[5]for fits A and B with fixed a s (m z ) = 0.118. 
As expected from the above discussion, curves for g 2 o = and g 2 o = g\% W are quite 
close to each other. Quark distributions are hardly affected by mass corrections. 
The polarized gluon density Ag is determined, in this procedure, only through the 
effect of scaling violations. It is therefore less constrained than quark distributions. 
However, the modifications of the gluon distribution induced by mass corrections, 
although sizeable, are considerably less important than the uncertainty originated by 
the choice of different functional forms for the input distribution (see ref. |[14|| ). 




Figure 5: The gluon distribution at Q 2 = 10 GeV 2 in fits A and B, with and 
without target mass corrections. 
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4. Conclusions 

We have computed target mass corrections to nucleoli structure functions in po- 
larized deep-inelastic scattering. Our results are consistent with those obtained in 
ref. [jTjJ using a different technique. 

Target mass corrections can in principle be important in the context of polarized 
deep-inelastic scattering, because part of the data are taken at relatively low values 
of Q 2 . For this reason, we have performed an analysis of all available polarized 
deep-inelastic scattering data in the framework of perturbative QCD, including the 
contribution of target mass corrections up to terms of order m 2 /Q 2 . We have proved 
that this approximation is reliable in the kinematical range of presently available data. 
We have found that the effect of mass corrections is generally small; for example the 
value of the axial coupling is enlarged by approximately 0.03 when target mass 
corrections are included. The strong coupling constant a s (mz) receives corrections of 
approximately 0.004. Both deviations are compatible with higher twist uncertainties 
estimated in previous works. 

Quark distributions are practically unchanged by the introduction of mass terms. 
The polarized gluon distributions, which is only determined through scaling viola- 
tions, is comparatively more affected by mass corrections, but also in this case the 
values of the fitted parameters do not deviate from those obtained for m = by more 
than one standard deviation. 
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